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Abstract 

Brownian Dynamics simulations are canied out to understand the effect of temperature 
and dielectric constant of the medium on microphase separation of charged-neutral diblock 
copolymer systems. For different dielectric media, we focus on the effect of temperature on 
the morphology and dynamics of model charged diblock copolymers. In this study we examine 
in detail a system of partially charged block copolymer consisting of 75% neutral blocks and 
25% of charged blocks with 50% degree of ionization. Our investigations show that due to the 
presence of strong electrostatic interactions between the chai^ged block and counterions, the 
block copolymer morphologies are rather different than their neutral counterpart at low dielec- 
tric constant, however at high dielectric constant the neutral diblock behaviors are observed. 
This article highlights the effect of dielectric constant of two different media on different ther- 
modynamic and dynamic quantities. At low dielectric, the morphologies are a direct outcome 
of the ion-counterion multiplet formation. At high dielectric, these charged diblocks behav- 
ior resembles that of neutral and weakly charged polymers with sustainable long-range order. 
Similar behavior has been observed in chain swelling, albeit with small changes in swelling 
ratio for lai^ge change in polarity of the medium. The results of our simulations agree with 
recent experimental results and are consistent with recent theoretical predictions of counterion 
adsorption on flexible polyelectrolytes. 

Introduction 

Microphase separation in neutral-neutral diblock copolymer melts and solutions has been stud- 
ied extensively during the last two decades.-"- However, self-assembly of amphiphilic macro- 
molecules still eludes a clear understanding. Despite the importance of these molecules for a 
number of technological and biological applications, a large parameter space affecting the self- 
assembly process poses serious problems for the scientific community. 



In particular, an amphiphilic diblock copolymer containing a charge block attached to a neu- 
tral block is of great importance to a large classes of nanotechnology applications -"— and drug- 
delivery systems.—"— Pioneer experiments on the self-assembly of amphiphilic copolymers in 
solutions by Eisenberg and his coworkers— have revealed that different morphologies can be ob- 
tained by varying the composition of polystyrene (PS)-Zj-poly(acrylic acid) (PAA). In addition to 
the three well-known classical morphologies, namely, spheres, cylinders or rods and lamellae, two 
more morphologies were discovered. These two morphologies included vesicles in aqueous solu- 
tion and simple reverse micelle-like aggregates in organic solvents. 

Origin of these morphologies (or micelles) in the context of amphiphilic copolymer solutions 
has been investigated extensively .-i^'^'^ General consensus is that the micellization/aggregation 
of an amphiphile results from a balance between three different contributions to the free energy: 
chain stretching in the core, the interfacial energy and the repulsion/attraction between the coro- 
nal chains.— However, the presence of charges in the case of charged-neutral diblock copolymers 
complicates this simple picture of the balance of forces and sometimes, leads to non-trivial counter- 
intuitive morphologies. For example, for a particular set of parameters, stacking disk-shape struc- 
tures have been observed experimentally .-^^-^^^ As an another example, in a system consisting 
poly(acrylic acid-Z7-styrene) (PAA-Z7-PS) diblock copolymer, disk-shape one dimensional supra- 
assembly have been formed in a controlled manner.—"^ 

To understand the structures and the mechanisms by which amphiphilic diblock copolymers 
self-assemble into different morphologies, a comprehensive theoretical/simulation study is neces- 
sary. However, theoretical work on the ordered morphologies of amphiphilic block copolymers is 
relatively scarce.—"— The few theoretical studies of diblock amphiphiles and triblock copolymers 
that have been conducted to understand the formation of self-assembled structures includes molec- 
ular dynamics (MD),-^^-^ Monte Carlo (MC)^"-^ and self-consistent field theory (SCFT).-^^'^ 
Although most of these recent simulation studies address self-assembly in solution,-^^^ recently 
there have been efforts to approach the problem using dissipative particle dynamics (DPD) to un- 
derstand the mesoscopic self-assembly of amphiphilic block copolymers in polymer melts.—"— In 



solution, Guo et.al.^^^ performed DPD simulation of paclitaxel loaded poly(ethylene oxide)-^?- 
poly(L-lactide) PEO-Z7-PLIA in water and dimethylformamide (DMF) to observe self-assembled 
structures of bicontinuous, lamella, rod and spherical micelles and provided a complete phase be- 
havior of the same system. Quite recently, it has been shown by Kriksin et.ah— that the local 
chemical structure of monomeric units can influence the global self-assembled morphology of the 
amphiphilic melts. It is interesting to note that so much effort has been invested to understand 
block copolymer assemblies, however, the explicit inclusion of charged blocks is by and large 
missing in most of the simulations of amphiphilic block copolymers. So far, the self-assembly 
of charged-neutral diblock copolymers has been studied using analytical techniques such as the 
random phase approximation-^"^ (RPA), the SCFT^^ and the MC.-^ In MC studies, ^^ only the 
thermodynamics of the system have been investigated due to the equilibrium limitations of the 
Monte Carlo technique which may do justice to understand the morphologies but lacks a thorough 
investigation that include dynamics. Only recently, experiments— have been performed to under- 
stand the microphase separation in the charged-neutral diblock copolymer melts. Theoretically, it 
has been predicted that the entropic cost of confining the counterions to charged domains as a result 
of microphase separation stabilizes the disordered phase against ordered ones for a large parameter 
range. Furthermore, the morphology diagram is highly asymmetric and the regime of stability of 
the lamellar phase gets enhanced. We must point out here that in the field theoretical models,—"^ 
possibility of counterion adsorption—"^ with the lowering of the temperature was not considered 
and hence, the models ignore the effect of ion pair formation and their implications on the struc- 
ture. This, in turn, limits the applicability of the field theoretical models without ion-pair formation 
to high temperature regime close to the disorder-order transition. 

In this article, we investigate the dynamics and self-assembly of charged diblock copolymers 
at nanoscale. The diblock copolymers consist of a charged block and an uncharged block. We are 
interested in the morphology and dynamics of charged block copolymers near the melt monomer 
density. Simulation of more realistic conditions requires an 'explicit' presence of the charge- 
counterion interactions. In order to allow for an efficient simulation but account for explicit 



Coulomb interactions between the charges we employ a Brownian dynamics (BD) technique with 
Kremer-Grest bead spring model polymer to understand the structural as well as dynamical prop- 
erties of the system and investigate the parameter space to identify the most important physics 
issues. The interplay between the entropy of the system and electrostatic energy between the 
charged monomers causes the formation of specific morphologies. We have observed that the di- 
blocks of 50% charge states on the charged block of an amphiphile forms structures that are quite 
different from their neutral counterpart. The chain swelling and other thermodynamic properties 
show the system does not follow the normal polyelectrolyte behaviors either. Similar results have 
been observed in recent experimental investigation s i ^i^i i^ for sulfonated polystyrene (sPS) and 
fluorinated polyisoprene (fPI) diblock copolymers. These new classes of charged block copoly- 
mers show promising new directions to understand and thereby fabricate novel functional materials 
which can be applicable to rather diverse applications from drug delivery to molecular electronics. 
The paper is organized as follows: Next section describes the simulation methodologies taking into 
account both Lennard- Jones (LJ) and Coulomb interactions. The relevancy of different parameters 
used in this simulation is explained in this section too. The results are discussed in third section 
and compared with existing simulation and experimental results. In sectiori, we conclude with a 
short description of the major findings and their importance to novel material designing. 

Simulation Method 

We used Brownian Dynamics (BD) simulations to examine the morphology and dynamics of 
charged block copolymer chains in a melt of density, p = Q.l /o^. The simulations are performed 
for two different system of sizes V = 16 x 16 x 16a^ and 2V = 20.16 x 20.16 x 20.16a^. For 
both the systems chain length of A'^ = 64 is used. Each chain contains 75% uncharged block and 
25% charged block. The charged block has 8 charges of charge +q on it, giving a 50% degree of 
ionization on the backbone. The charges interact via Coulomb forces with equal number of counte- 
rions randomly dispersed in the system. The initial configuration of the model system is randomly 



generated at the melt density. All the monomers in the system have mass m,- and Lennard- Jones di- 
ameter, o. Polymer chains are modeled following the Kremer-Grest bead spring polymer model- 
in which bonded beads are connected by finitely extensible nonlinear elastic (FENE) springs rep- 
resented by, 
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where Rq = 1.5a is a finite extensibility and the spring constant, k = 37.5£/a^, o being the 
monomer diameter. The FENE potential in combination with the (excluded volume) repulsive in- 
teraction creates a potential well for the flexible bonds that maintain the topology of the molecules. 
The energetic interaction between the beads is modeled by a truncated and shifted Lennard- Jones 
potential. The repulsive part of the potential is given by. 
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And the attractive part of the potential is given by. 
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where and rjj is the distance between two monomers and £r and £a are the repulsive and 
attractive energy parameters respectively for two different interactions described below. Each 
monomer of the system interacts via a short-range repulsive potential whose interaction strength, 
£ ^ £r = 1.0. The short-range repulsive LJ potential is shifted and truncated with a cut-off dis- 
tance, rf: < 2^'^a. In addition to this repulsive interactions, blocks A-A and B-B are attractive to 
each other with interaction strengths, e^ = 2.0 and £bb = 4.0 respectively. The above choice de- 
rives from the fact that the natural tendency of the different blocks of a BCP is to avoid each other 
because of the presence of dispersive intermolecular forces which often results in similar blocks 



more attractive to each other than the dissimilar monomers.— The cross interactions for A-B with 
strength, £ab, are then obtained using Lorentz-Berthelot mixing rules. The repulsive cutoff is used 
in conjunction with the attractive cutoffs for the AA and BE interactions with a cut-off distance, 
^attractive ^ 2.5a. As the focus of this paper is to observe the effect of electrostatics and entropy on 
the self-assembly and dynamics of the chains, we vary the temperature of the system along with 
the electrostatic parameters. 

The long-range interactions between the charges on the chain and the counterions are modeled 
using explicit Coulomb potential, 

where D is the dielectric constant of the medium. D may not be constant throughout the sys- 
tem if dielectric mismatch is considered, however, for the purpose of this paper D is assumed 
to be spatially constant. Long range Coulomb interactions are accounted for through the Ewald 
sum.— Temperature is the energy scale parameter that is varied for different sets of simulations. 
We introduce a second energy scale parameter which is the ratio of Coulomb and Lennard- Jones 
interactions: ^b = Q^ /{Do£r). For real experimental systems, ^b ~1-100: we use two values of 
t,B = 2 and 10. The parameter E,b is proportional to the Bjerrum length, /g and is a constant. 
The dynamics of the monomers are governed by the classical Newton-Langevin equation, 

m§ = -^U,^r§^W,(,) (5) 

where Ui is the net potential energy experienced by particle / and m, is its mass. F is the fric- 
tion coefficient between the chain monomer and background solvent. Wi{t) represents a Gaussian 
'white noise' with zero mean acting on each particle.—'^ The last two terms couple the system to 
a heat-bath where the 'friction term' acts as a heat sink and the 'noise term' acts as a heat source. 
The first advantage of this scheme is that the natural MD integration time-steps are larger, thereby 
permitting simulation on longer time scales. A second advantage comes from the fact that on this 
time-scale, only the mean effect of the stochastic forces acting on the system needs to be con- 



sidered, leading to the first order temperature relaxation law which in turn reduces the need of an 
external thermostat. The dimensionless units are defined as follows, f^ = t j \fm^^JeR, p* = po^, 
T* = keT/e, U* = U /keT and r* = r/o. 

Results and Discussion 

In earlier studies ,-22j^ we have demonstrated that controllable morphologies can be achieved by 
changing charge states of the charged diblock copolymer or by changing the dielectric constant of 
the solvent. In an attempt to establish a connection between thermodynamics to the dynamics of 
the system, we have performed a detailed molecular dynamics simulation of charged polymer and 
investigate the dynamics of the BCP and counterions for different thermodynamic parameters. To 
understand the system size effect we studied two systems, V and 2V described earlier. We also have 
varied the Coulomb energy parameter, i^g. The interaction parameter, ^b is inversely proportional 
to the dielectric constant, therefore t,B = 10.0 and 2.0 are used to represent low and high dielectric 
constant systems respectively. The system size effect is studied for i^b = 10 cases only. 

[Figure 1 about here.] 



[figure] [1][] 1 shows the morphologies of the block copolymer melt for 2V system at ^b = 10. 
In these figures, the conventional diblock morphology as expected from the phase behavior of 
neutral block copolymers cannot be observed. From the morphology diagram for neutral diblock 
copolymers ,-i^ it is well established that the neutral counterpart of this system of 75-25 diblock 
would have exhibited hexagonally packed cylindrical morphology. In case of neutral diblock, the 
minority components form the hex structures and the majority blocks form the bulk. The presence 
of charge sites causes electrostatic interactions to dominate over the energetic repulsion effects be- 
tween blocks that induces microphase separation in the neutral BCP. The strong electrostatic effect 
in the charged BCP changes the interaction energies such that the hexagonal structures cannot be 
observed. Lowering the temperature to T* = 0.3, gives rise to 'inverse' morphology where the mi- 



nority blocks form the structures and the majority components form the matrix ( [figure] [1][]1 ^b) 



At even lower temperature, T* = 0.05, the cylindrical nature of the 'inverse' morphologies can- 
not be observed as can be seen in [figure] [1][]1| ^ a). At the higher side of the temperatures, the 



charged blocks start to form well defined structures ([figure] [!][]! c)) that eventually settles down 



to entropy dominated miscible phase at very high temperature, T* = 1.0 ( [figure][l][]l| ^d)). The 



structures seen in [figure] [1][] 1 have also been observed in recent experiments.— The strong elec- 
trostatic interaction causes percolation of the charges that forms a network of charged blocks at low 
temperatures, giving rise to these unconventional 'inverse' morphologies. Furthermore, note that 
the counterions are confined in the charged domains for all the temperatures corresponding to the 
ordered morphologies, in agreement with the field theoretical calculations on similar systems.— 
The lower the entropy of the system, stronger the effect of the electrostatics, giving freedom to 
the uncharged block to self-assemble to novel morphologies that would not have otherwise been 
possible. 

[Figure 2 about here.] 



Snapshots for high dielectric constant at ^b = 2.0 are shown in [figure] [2] []2 Interestingly 
these snapshots show similar behavior as neutral 75-25 diblock. At lower temperature end, in 
[figure][2][]2| ^a) and (b), the charged block (minority component) form the micro structures and 



the majority neutral blocks form the bulk. By increasing the temperatures the structures break 
down as can be seen in [figure] [2] []2| ^c) and (d) for T* =0.8 and T* = l.O respectively. At high 



dielectric constant (low ^b), the agglomeration of charges is more prevalent compared to low D 
(high £,b). In a recent experiment, "^^ it has been shown that the long-range structures of charged 
diblock copolymers (sPS-fPI) breaks down by adding a trace amount of water in these system. It 
is well established for heterogeneous materials that an increase in dielectric constant results in a 
decrease in percolation and vice versa.— It is interesting to observe the same physics dominat- 
ing the microphase separation of macromolecules, in this case, the charged diblock copolymers. 
In [figure][2][]2| ^a) and (b), increase in dielectric constant resulted in reduced percolation thereby 



allowing the charges (minority charged block) to agglomerate to form the structures whereas the 
majority block forms the bulk. A detail discussion on agglomeration of charges will be presented 



in connection with radial distribution function and fraction of free counterions later in this section. 
These results are in agreement with experiments^ and the SCFT calculations— carried out for 
weakly charged diblock copolymers. The addition of water increases the dielectric constant of 
the system which in turn lowers the electrostatic interaction strengths by a factor of IJD. This 
makes the percolation (responsible for 'inverse' morphology) of charges more unstable, resulting 
in charge agglomeration (responsible for neutral BCP morphology) instead. Decreasing the elec- 
trostatic strength brings the system from an electrostatic energy dominated regime to an entropy 
dominated regime, where ion-pair formation is weaker compared to a percolated state. In the en- 
tropy dominated regime, the microphase separation is hindered by the entropic cost of partitioning 
of counterions in the charged domains as predicted by the SCFT calculations.— This also explains 



the breakdown/disappearance of 'inverse' morphologies in [[figure] [2] [] 2| obtained at low t,B- In 
principle, one should obtain the 'inverse' morphologies even for high D {^g = 2) medium by low- 
ering the temperature further. However, it is extremely time consuming and hard to equilibrate the 
system at such low temperatures (with explicit Coulomb interactions). 

[Figure 3 about here.] 

In order to corroborate the idea of ion-pair formation and its implications on the morphologies. 



we show the counterion-counterion radial distribution function (RDF), g{r), in [figure][3][]3 For a 
complete understanding of ion-pair formation, a discussion of the ion-counterion RDF is extremely 
important. Due to the lack of data of fixed charges (co-ion), the ion-counterion RDF cannot be pre- 
sented in this paper. Recent work on charged polymers— shows that the counterions always pair 
with ions. Ion-counterion pairing is also necessary from the electrostatics point of view, although 



this pair formation does not guarantee local electroneutrality. In [figure] [3] []3 'a) and (b) the RDF 
for ^5 = 10 and ^b — ^ are shown respectively for the system size, V. We confirmed by plotting 
the RDF for 2V system that there are no apparent system size effects which might cause substantial 
error in the results presented in this paper. The system size effects for the simulation model can 
be neglected with confidence. One of the prominent features of these radial distribution functions 
can be seen by comparing the RDF for different temperatures. The first peak height increases by 
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decreasing the temperature which represents increased agglomeration of counterions. The RDF 
shows a distinguishable decrease in agglomeration peak for (^6 = 2 compared to <^b = 10. This rep- 
resents an enhanced miscibility shown in [figure] [2] []2 In [figure] [3] []3| ^a) for low D, structures 



can be observed at temperatures as high as T* = 1.0 which is prominent in the snapshots shown 
in 



[figure] [1][]1 1 For t,B = 2 system, an order-to-disorder like transition can clearly be observed 



from the sudden decrease in the g{r) peak in |[figure][3][]3| ^b). The transition temperature is hard 



to locate exactly, these plots can give a qualitative understanding of transition that is associated 
with charged diblock copolymers. In polymer melt, it had been shown— that the precise transition 
temperature was very hard to locate. 

For the model of a continuously and uniformly charged line, the radial distribution function 
implicitly has two peaks, the first peak corresponds to the condensed counterions and a second 
peak corresponds to the Debye-Huckel peak for the uncondensed counterions. The first peak is 
obviously an idealization of the real structure. In [figure] [3] []3| ^a) and (b), the appearance of two 



peaks in RDF at low temperatures represents the counterion adsorption behavior akin to uniformly 
charged polyelectrolytes. However, this is not quite true because of the presence of a longer un- 
charged block, electrostatics and entropy both play a major role in deciding the morphology of 
the charged BCR The increase in T* breaks the morphologies at high D ( [figure] [2] []2[ ) which 



can also be seen in [figure] [3] []3 'b) where the first peak height is reduced. In [figure] [3] [] 3 a) 



the peak heights are higher and represent stronger agglomeration compared to the high dielectric 



constant case, [figure] [3] []3 'b) which has lower peak heights representing a much weaker agglom- 



eration, in this case the entropic effect is much more pronounced in forming the morphologies. 



In [figure] [3] []3 'a), the second peak becomes dominant which may result in longer inter-atomic 



separations in a low dielectric constant medium. The first peak in [figure] [3] []3| ^a) begins dis 



appearing and the majority of counterions show agglomeration at an uncondensed phase (second 
peak) akin to the polyelectrolyte behavior. Also at high T* the thicker width of ^(r) is represen- 
tative of loosely bound counterions similar to polyelectrolyte behavior.— In|[figure][3][]3|^b), the 



first peaks shift to shorter inter-atomic separations with the second peak position unchanged. For 
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high dielectric medium, ^b = 2, the first peak positions shift due to tighter packing. The close 



packing of counterions seen in [figure] [3] []3 can electrostatically be possible if the counterions 
form charge-counterion pairs. Therefore, the neutral diblock behavior of the morphologies as ob- 
served in large D ( [[figure] [2] []2 ) are due to the the tightly packed charge-counterion pair at each 



charge site of the polymer that causes the charge sites to resemble as neutral sites. On the other 
hand, the presence of the second peak indicates counterion-counterion correlations that extends 
throughout the charged block of the chain which effectively shows apparent neutralization of the 



charged block. This relates quite well with the snapshots shown in [figure][l][] 1 and [figure][2][]2 
From these morphological analysis and the structural evolutions for different dielectric constants, 
it can easily be concluded that at low D the charged block exhibit polyelectrolyte behavior and 
at high D it behaves like a neutral diblock. It should be noted that both the ^b studied here have 
explicit point charges on the minority block of the chain. 

[Figure 4 about here.] 

To augment the understanding of charge agglomeration, cluster size distribution is investigated. 
A cluster is defined to be the total number of nonbonded charges associated with 'a' charge within 
a cutoff radius rcut- The cut-off radius, rcut is the location of the first minimum of RDF in the max- 
imally clustered state. The total number, Nc, of non-interacting charge sites within rcut is defined 
as the cluster size. The probability of occurrence, P{Nc), of a cluster of size Nc is computed. The 
probability of occurrence of cluster size A^c = 1, ^(1) is taken as the fraction of free counterions. In 
[figure] [4] []4[ fraction of free counterions are plotted versus temperature for two different ^b- The 



red squares represnt i^g = 10 and black circles represent i^s = 2 respectively. Small P(l) reflects 
lower free counterions, threfore higher agglomeration. For lower temperatures, very few counte- 
rions are free. Stronger agglomeration can be observed in snapshots at lower temperatures too. 
There is a sharp increase of /'(I) at around T* =0.1 which may be associated with stronger ag- 
glomeration of charges. For ^b = 10, the red sqaures show more clustering (less free counterions) 



compared to ^b = 2. From the snapshots (|[figure][l][]l|and|[figure][2][]2[), it has been observed 



that low dielectric constant (^b =10) shows more percolation of charges thereby forming reverse 
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morphologies. These resuks are consistent with [figure] [4] []4 in which it can be observed the 
^B = 10 forming more clusters (less free counterions). From these plots though the formation of 
charge multiplets cannot be guarateed as that requires an explicit representation of the fixed charge 
data which is absent due to unavailability of fixed charge postions. However, agglomeration of 
counterions cannot be justfied, from a electrostatics point of view, without the presence of oppo- 
site charges. For charged polymer systems, it has been observed earlier— that multiplets form in 
charge agglomeration. 

[Figure 5 about here.] 

From these results presented so far, it is evident that thermodynamics of the charged-neutral 
diblock copolymers depend strongly on the electrostatic interaction strengths. To understand the 
thermodynamics better we have computed the specific heat capacity (Cy) from the fluctuation- 
dissipation formula, Cy = {< E^ > — < E >^)/T*^ where E is the total energy of the system and 
'<>' denotes ensemble average. [figure][5][]5| ^a) shows the Cy as a function of temperature for 



different ^b- For ^g = 10 (red squares), a peak in Cy is observed at around T* = 0.1 at low dielectric 
constant. The peak is absent for i^b = 2 (black circles) at higher dielectric constant. Instead there 
is sudden jump in Cy for high dielectric constant {^b = 2)at around T* = 0.1. This suggests that 
the charged blocks undergo self-assemble (order-to-disorder transition) at this temperature. At 
around this temperature, the entropy of the system decreases thereby causing an abrupt change in 
Cy. Typically, the appearance of a 'peak' (red squares) in Cy is often referred to as the "glass 
transition temperature". However, the peak observed here for ^b = 10 is not associated with a 
glass transition as the system cannot vitrify at this temperature and density. At i^g = 10, the self- 
assembly is driven by electrostatic interactions rather than by entropic interactions. Since the Cy 
peak occurs at low dielectric constant, it is evidently related to the formation of ion-counterion 
multiplets. ~ The multiplet formation starts occurring much earlier than this transition temperature 



as is evident from | [figure] [ 1 ] [] 1] This phenomena cannot be observed at i^g = 2 as the charges are 
dissociated at high dielectric constant. At i^g = 2 the change in Cy comes from the conformational 
energy of the chain rather than the multiplet formation as is evident from the abrupt increase in 
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Cy at low T* . Furthermore, the absence of an abrupt increase in Cy for ^b = 10 (low dielectric 
constant medium) can be attributed to the higher order nature of the transition that may be caused 
by strongly percolating charged blocks. We have not investigated the exact nature of this transition 
and leave this finer point for future work. The structures also break down at higher dielectric 
constant. The same phenomena has been observed in experiments^ where an addition of small 
trace amount of water (increasing dielectric constant) breaks down the long range order of the 
charge block copolymers. 

To understand the structural changes that generally affect the thermodynamics, we have com- 
puted the structural relaxation time of counterions and plotted in [figure] [5] []5| ^b). The first dy- 



namical quantity in this article, the structural relaxation time is obtained from self-intermediate 
scattering function defined by, S{k*,f') where k* is a wave vector, k = 2n/l corresponds to the 
first peak distance (/) of the static structure factor, S{k*). The intermediate scattering function, 
S{k*,t* follows an exponential decay. The structural relaxation time, T* is represented by the value 
off* where S{k*,t*) resumes l/e times initial amplitude. The red squares and blue triangles are 
shown for counterions and chain center-of-mass (CM) respectively for ^b = 10^ i-e-. at low dielec- 
tric constant and the black circles are shown for counterions at i^g = 2. For ^b = 2, the chain CM 
relaxation time is not shown as the major focus of this investigation is the effect of counterion re- 
laxation. At low dielectric constant, formation of charge multiplets does not allow the counterions 
relax thereby a slow relaxation can be observed (red squares). These plots show an abrupt change 
in T* at low temperature which is consistent with the Cy plot shown in [figure] [5] []5| ^a). For the low 



dielectric constant {^b = 10), the counterions cross the chain CM at around T* =0.1, exactly the 
temperature at which the peak in Cy has been observed. These results support the theory that the 
counterions undergo a structural transition at low temperatures due to multiplet formation. Below 
this temperature the counterion relaxation time increases strongly, consistent with the localization 
of charged blocks into multiplet clusters. A similar behavior has been observed in other weakly 
charged systems where the transition occurs due to the formation of charge multiplets.— 

[Figure 6 about here.] 
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To understand the swelling behavior of the charged diblock copolymer for different dielectric 
constants, we compare the radius of gyration of the charged polymer (Rg) with that of an ideal 



chain, Rgo = \/Nb'^/6, where N is the degree of polymerization and b is the monomer length. 
The swelling ratio, 0^ = Rg/Rgo versus temperature is shown in [figure] [6] [] 6 The ratio, 0^ 



is similar around the transition temperature for both E,b = 2 (black circle) and 10 (red squares). 
However, beyond the transition temperature, 0^ is slightly higher for (^5 = 10 than ^b = 2. This 
can be explained as follows. As the dielectric constant increases (i^g = 2), swelling capability of 
the chain decreases which can normally happen in neutral polymers only. Therefore, the swelling 
phenomena of the charged BCP system for i^b = 2 may be attributed to the swelling of neutral 
chains. On the other hand, for polyelectrolytes, increase in dielectric constant increases swelling 
capability. Therefore for i^g = 10 (red squares), increase in swelling phenomena can be attributed 
to the polyelectrolyte behavior of the charged BCP. In these data comparatively small swelling 
of the chain has been observed, in some cases though, the swelling has been observed to be as 
much as 100 times their volume for particular classes of polyelectrolytes in low-polarity (low- 
dielectric constant) solvents.— This suggests that, although there are dissociation of counterions 
to break the multiplets with the increase in dielectric constant, the entropy due to the dissociation 
of counterions is relatively small to contribute substantially in swelling of the charged BCP. The 
system under consideration has only 8 charges (and 8 counterions) on the first 16 monomers of 
a chain of length 64. These small number of charges result in statistically smaller contribution 
to the entropy of dissociation, although the Coulombic interactions significantly contribute to the 



formation and breaking of charge agglomeration as shown in [figure][l][]l and [figure] [2] []2 

[Figure 7 about here.] 
[Figure 8 about here.] 



[figure] [7] []7 and [figure] [8] [] 8 show the mean-square-displacement (MSD) of the counterions 



and chain CM for two different t,B. For i^g = 10 ([[figure] [7] []7[), only the V system is shown due 



to lack of data to calculate correlation in the 2V system, but this should not be a major concern 
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as there is hardly any system size effects present in these calculations. In [figure] [7] []7| ^a) and 



[figure] [8] [] 8 counterion motion is observed to be much faster than the chain CM motion (shown 



m 



[figure] [7] []7^b) and [figure][8][]8;b)). At low temperatures, both counterions and chain CM 
show sharp decrease in their respective MSDs. Although counterions show faster motion due to 
their smaller size compared to the chain CM, the drop in counterion MSD is more pronounced than 
the chain CM at T* = 0.05 and 0.1 respectively. This may be due to the stronger agglomeration of 
counterions that causes the counterions to be entropically unfavorable to move. For low dielectric 



constant, i^g = 10 (|[figure][7][]7|^a) the low temperature counterion motions become even more 



sluggish. This is consistent with the earlier observation of counterions forming charge multiplets 
at low temperatures for high ^b- For ^b = 2, the counterions dissociate due to the presence of more 
polar medium (high dielectric constant) therefore the motion is enhanced, although slow compared 
to the chain CM. The sudden jump in MSD for both the chain CM and counterions for both (^5 = 2 
and 10 suggests a transition from the non- agglomerated state to charge agglomerated state. 

Conclusion 

In summary, we have simulated a system of charged block copolymers in the presence of explicit 
counterions. The partially charged block copolymers exhibit 'inverse' morphologies that have 
been observed in recent experiments.— From our results, it is evident that the design of these 
novel functional materials depend sensitively upon the dielectric constant of the medium, the de- 
gree of ionization and the temperature of the system. Corroborative results have been obtained 
experimentally— that addition of trace amount of water (increasing dielectric constant) changes 
the morphology of the diblock copolymer dramatically. 

Furthermore, our results highlight the existence of two regimes in the self-assembly of charged- 
neutral diblock copolymers. In one regime, the self-assembly is governed by electrostatic energy 
and in the other regime, it is dictated by the counterion entropy affecting the morphologies and 
other thermodynamic properties. It is interesting to note that introducing charges on a diblock 
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copolymer not only changes the morphologies, it acts unfavorably for both neutral polymers and 
polyelectrolytes depending on the dielectric of the medium. In a low dielectric medium, the 
self-assembly in charged-neutral diblock is electrostatically dominated and leads to non-trivial, 
counter-intuitive 'inverse' morphologies where the minority component forms the matrix and ma- 
jority forms the structures. For this system, electrostatics drive the formation of ionic multiplets 
which increases the relaxation time of the counterions and hence affect the dynamics of microphase 



separation ( [figure] [1][] 1 ). For high dielectric medium, the chains form neutral BCP morphologies 



as can be seen in [figure] [2] [] 2 An increase in dielectric constant shifts the system towards coun- 
terion entropy dominated regime. The counterion dominated regime also can be observed with 
increase in temperature. In this regime, disordered phase gets stabilized against ordered structures. 
These results exhibit how different the 75-25 charged BCPs are from the neutral diblocks coun- 



terparts. The swelling behavior shown in [figure] [6] [] 6 shows how different these charged BCPs 



are from their polyelectrolyte counterparts. The swelling of the chains, although very small, shows 
behavior which is opposite to polyelectrolytes swelling behavior at low D. The formation of charge 
multiplets in a cluster and the dissociation of counterions with decreasing dielectric constant, can 



be understood by observing the dynamics of counterions ( [figure] [7] []7 and [figure] [8] [] 8 1. These 
results are in agreement with the SCFT calculations^ for weakly charged systems, showing the 
stabilization of disorder phase against ordered morphologies to avoid the entropic cost of partition- 
ing of counterions in the charged domains. It should be noted here that a discussion on the slope 
of the MSD (diffusion coefficients) is missing. A detailed discussion of the diffusion coefficients 
of the chain CM and counterions may indeed be important to understand the effect of transport 
coefficients on the morphological changes for different T* and D, however, diffusion in polymers 
is a subject that deserves a complete separate study. 

The present structural and dynamical analysis of charged diblock polymers should be of im- 
portance for understanding the fundamental mechanisms and for designing novel materials for 
drug delivery I3ii4i2i,24 ^ ^^^ nanotechnology— "^ applications. Additional information of differ- 
ent charged states and block sizes of the charged block will further enhance the understanding and 
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in particular provide information for counterion transport in charged polymer films. 
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List of Figures 



Snapshots at the end of the simulation for different temperatures at i^g = 10.0 for 
2V system, (a) T* = 0.05, (b) T* = 0.3, (c) T* = 0.8 and (d) T* = 1.0. The maroon 
represents uncharged block monomers. Black and yellow represent charged block 
monomers, black being the neutral monomeric unit and yellow being the charged 
sites of the charged block. The green dots are the counterions. The central sim- 
ulation cell is repeated in all the three directions to preserve the continuity of the 

chains images due to the periodic boundary conditions |25 

Snapshots at the end of the simulation for different temperatures at ^b = 2.0. (a) 
T* = 0.05, (b) T* = 0.1, (c) T* = 0.3 and (d) T* = 0.8. The maroon represents 
uncharged block monomers. Black and yellow represent charged block monomers, 
black being the neutral monomeric unit and yellow being the charged sites of the 
charged block. The green dots are the counterions. The central simulation cell is 
repeated in all the three directions to preserve the continuity of the chains images 

due to the periodic boundary conditions |26 

Radial distribution function g{r) of counterions for all the temperatures. The two 
different systems shown here are: (a) ^b = 10 at V system size and (b) i^s = 2 for 
V system size. For ^b — 2, more temperature data have been obtained as shown in 

the legend |27 

Fraction of free counterions as a function of temperature for t,B = 2 (black circles) 
and t,B = 10 (red squares) systems. The lines in these plots are guide to the eyes. . |28 
Specific heat, Cy for all the temperatures T* = 0.05, 0.1, 0.3, 0.5, 0.8 and 1.0. (b) 
Counterion structural relaxation time for i^^ = 10 (red squares) and i^b = 2 (black 
circle). The blue triangles represent polymer center-of-mass structural relaxation 
time for 2V system at ^b = 10. At around T* = 0.3 the counterion T* becomes 
stronger than the chain CM T* l29 
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Ratio of radius -of -gyration to the radius-of-gyration of an ideal chain of same chain 
length for t,B = 2 (black circles) and ^5=10 (red squares) systems. For t,B = 10, 
system size 2V is used, ^b = 10 shows more chain swelling compared to ^b = 2. 
The lines are not data, these are merely guide to the eyes for the respective swelling 

parameters |30 

Mean square displacement at ^b = 10.0 for (a) counterions and (b) chain center- 
of-mass. Black line represent T* = 0.05, red line 0.1, green line 0.3, blue line 0.5, 
magenta line 0.8 and cyan line 1.0 respectively. MSD is shown for V system due 

to lack of data for 2 x V system |31 

Mean square displacement at i^g = 2.0 for (a) counterions and (b) chain center-of- 
mass. Respective temperatures for the lines are shown in the legends \32 
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(a) T* = 0.05 



(b) T* = 0.3 





(c) r*=o.; 



(d) r* = i.o 



Figure 1: Snapshots at the end of the simulation for different temperatures at ^b — 10.0 for 2V 
system, (a) T* = 0.05, (b) T* = 0.3, (c) T* = 0.8 and (d) T* = 1.0. The maroon represents 
uncharged block monomers. Black and yellow represent charged block monomers, black being the 
neutral monomeric unit and yellow being the charged sites of the charged block. The green dots 
are the counterions. The central simulation cell is repeated in all the three directions to preserve 
the continuity of the chains images due to the periodic boundary conditions. 
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(a) T* = 0.05 



(b) r*=o.i 





(c) r*=0.3 



(d) T* = 0. 



Figure 2: Snapshots at the end of the simulation for different temperatures at ^b = 2.0. (a) 
T* = 0.05, (b) T* = 0.1, (c) T* = 0.3 and (d) T* = 0.8. The maroon represents uncharged 
block monomers. Black and yellow represent charged block monomers, black being the neu- 
tral monomeric unit and yellow being the charged sites of the charged block. The green dots are 
the counterions. The central simulation cell is repeated in all the three directions to preserve the 
continuity of the chains images due to the periodic boundary conditions. 
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(b) ^B = 2 

Figure 3: Radial distribution function g{r) of counterions for all the temperatures. The two dif- 
ferent systems shown here are: (a) ^b = 10 at V system size and (b) i^b = 2 for V system size. For 
(^B = 2, more temperature data have been obtained as shown in the legend 
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Figure 4: Fraction of free counterions as a function of temperature for i^s = 2 (black circles) and 
<^B = 10 (red squares) systems. The lines in these plots are guide to the eyes. 
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(b) Structural relaxation time 



Figure 5: Specific heat, Cy for all the temperatures T* = 0.05, 0.1, 0.3, 0.5, 0.8 and 1.0. (b) 
Counterion structural relaxation time for (^^ = 10 (red squares) and i^s = 2 (black circle). The blue 
triangles represent polymer center-of-mass structural relaxation time for 2V system at <^s = 10. At 
around T* = 0.3 the counterion T* becomes stronger than the chain CM T*. 
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Figure 6: Ratio of radius-of-gyration to the radius-of-gyration of an ideal chain of same chain 
length for t,B = 1 (black circles) and i^g = 10 (red squares) systems. For i^g = 10, system size 2V 
is used, (^g = 10 shows more chain swelling compared to t,B = 2. The lines are not data, these are 
merely guide to the eyes for the respective swelling parameters. 
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(b) Chain CM MSD 

Figure 7: Mean square displacement at i^g = 10.0 for (a) counterions and (b) chain center-of-mass. 
Black line represent T* = 0.05, red line 0. 1 , green line 0.3, blue line 0.5, magenta line 0.8 and cyan 
line 1 .0 respectively. MSD is shown for V system due to lack of data for 2 x V system. 



31 



10" 



10 



Q 2 

00 10 



S 10 



0.1 



r 

r 

r 


1 


*=^^ 


— T*=0.05 

— T*=0.1 

— T*=0.3 

— T*=0.5 

— T*=0.8 
T*=1.0 


^ 


r 


^^ 


^-^^^^^^^^^^ 


r 


^-^ 


^^^^^ 1 


■ 


-rTTTiTTi 


" 



10 100 

IntJr't*^ 



1000 



10^ 



10 



P 

00 



10- 



M 1 

O 



10" 



lO" 





1 1 1 








- T*=0.05 


^ - 


. 


- T*=0.1 


^/^^^ ■, 




- T*=0.3 


^^^^^^ 


■ 


- T*=0.5 


^<^^^^^^ - 


. 


— T*=0.8 


^-•ii-^^^^'^ ^-^'^ ~ 




T*=1.0 


^^j;^^^^^^^,^'^^^ - 


^ 


^^^^^ 1 


^^^^ 


.^^^ ■? 


c^^^^^^ 




1 





10 100 

log(t*) 

(b) Chain CM MSD 



1000 



10000 



Figure 8: Mean square displacement at h,s = 2.0 for (a) counterions and (b) chain center-of-mass. 
Respective temperatures for the lines are shown in the legends. 
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